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I. INTRODUCTION 

The detailed study of the properties of unstable nu- 
clei far from the stability line has been one of the main 
goals of nuclear physics in the last years. Recently, the 
development of radioactive ion beam (RIB) facilities in 
GSI (Germany) and in RIKEN (Japan) has opened a new 
field for such study, making possible the production of a 
variety of exotic nuclei which may have large neutron or 
proton excess. 

So far studies have largely dealt with the light nuclei 
which became accessible by the use of RIB produced 
in fragmentation reactions. Prior experiments have re- 
vealed a halo phenomenon (e.g. in Refs. 0, 0|) whose 
occurrence is due to both the small separation energy of 
the last few nucleons and their occupation on the orbits 
with low angular momentum. As well as halos, heavier 
systems may lead to formation of pronounced neutron 
skins [3] due to availability of far more neutrons than 
protons in these nuclei. 

Immediate determination of the neutron skin thick- 
ness usually involves the precise measurement of the root 
mean square (rms) radii of both charge and mass distri- 
butions. Electron- nucleus scattering has proven to be 
an excellent tool for the study of nuclear structure. In 
particular, it has accumulated much reliable information 
on charge density distributions of stable nuclei. There- 
fore, it is believed that the new facilities in GSI 0, Q 
and RIKEN @ will provide a good opportunity to study 
the charge density, and consequently the proton density 
distribution, of unstable nuclei by elastic electron scatter- 
ing. Unfortunately, a measurement of the neutron den- 
sity distributions to a precision and detail comparable to 
that of the proton one is hardly possible. The nuclear 
matter distribution in 6 He and 8 He has been determined 
recently at GSI by using small angle proton scattering 
in inverse kinematics at relativistic energy Q, and data 
has also been collected for n Li. It turned out that to 
get information on the neutron skin thickness one needs 



data obtained with probes having different sensitivities 
to the proton and neutron distributions. The methods 
for extracting the neutron skin thickness mostly include 
hadron scattering [3, Q , antiprotonic atoms , parity 
violating electron scatteri ng fill . [l2l [l3| , as well as giant 
dipole resonance method |l4[ and spin-dipole resonance 
method [HG1. 

On the theoretical side, calculations of nuclear charge 
and matter radii of exotic nuclei are usually made in 
the framework of mean-field approaches, namely Hartree- 
Fock (HF) method (see for example [13, [TH) or Hartree- 
Fock-Bogoliubov (HFB) method including pairing corre- 
lations [ij, S3, S3, S3, S3, Si- The latter predicts well 
the monotonic increase of the neutron skin thickness for 
different chains of isotopes up to the drip line [l9[ • Re- 
cently, the self-consistent relativistic mean-field (RMF) 
model has been widel y ap plied to both stable and unsta- 
ble nuclei (e.g., Refs. [H S3, S3, SI)- Also the relativis- 
tic Hartree-Bogoliubov approach has been employed to 
study the nuclear skin thickness in neutron/proton-rich 
sodium isotopes (23 |. Many calculations show that the 
RMF model can reproduce with good precision a number 
of ground-state nuclear properties including the charge 
radii [3fj |. The charge rms radii were successfully de- 
scribed very recently in Ref. , where a generator coor- 
dinate method (GCM) on top of Gogny HFB calculations 
was explored. 

Theoretical identification of skin and/or halo structure 
in neutron-rich weakly bound nuclei, however, is still a 
matter of discussion. In Ref. [3l| a definition of the neu- 
tron skin and its appearance were presented in terms of 
spherical HF calculations. The proposed criteria which 
deal with proton and neutron densities allowed one to 
predict neutron skins in nuclei far from the [3 stability 
line. It has been also shown in [3l[ that the formation 
of proton skin appears to be rather difficult. The Helm 
model [H, S3] has been applied in Ref. ,20?] to analyze 
neutron and proton skins, as well as halos, of even-even 
Ni, Sn, and Pb isotopes in terms of form factors. In Ref. 
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[20| three different definitions were proposed for neutron- 
proton radii differences. Among them, the one based on 
the Helm model was chosen as a measure of the skin. 
The latter was shown to have a smooth gradual depen- 
dence on the neutron excess and to be almost unaffected 
by shell effects. The Helm model was used very recently 
also by Bertulani [34| to investigate electron scattering 
from light unstable nuclei. 

Hitherto, the different definitions for the skin thickness 
mentioned above have been explored within different nu- 
clear structure models. We would like to emphasize that 
a comparison of skins extracted by using various defi- 
nitions is not very meaningful unless the same nuclear 
model is used and this has not been done in the past. 
Such an analysis of neutron skins within a given micro- 
scopic nuclear structure model could be very useful also 
in respect to demonstrate their expected spreading when 
different definitions of the nuclear skin are used. 

Another interesting question is to explore how the neu- 
tron skin emerges in the presence of deformation. The 
latter is defined by the non-spherical components of the 
proton and neutron density distributions. In particular, 
studies of deformed exotic nuclei and skins can be found 
in Refs. [H, Hfil ■ It is desirable to study the evolution of 
shape and skin formation, not only because deformation 
influences the nuclear rms radii, but also because of the 
possible skin anisotropics that may take place. 

In the present study, the properties of even-even Ni 
(4=48-78), Kr (A=70-I00), and Sn (4=100-136) iso- 
topes are described using the deformed self-consistent 
mean-field Skyrme HF+BCS method. We have used 
three parametrizations of the Skyrme force, namely SG2, 
Sk3 and SLy4, which were able to give an appropriate de- 
scription of bulk properties of spherical and deformed nu- 
clei in the past. As in our previous paper [22|], we choose 
some medium and heavy Ni, Kr, and Sn isotopes because 
many of these sets, which lie in the nuclear chart between 
the proton and neutron drip lines can be formed as ra- 
dioactive ions to perform scattering experiments. The 
main goal of this study is to clarify theoretically the emer- 
gence of the neutron and proton skins in neutron-rich and 
neutron-deficient isotopes, respectively, by testing differ- 
ent definitions for the skin thickness in the framework 
of the deformed Skyrme HF+BCS model. Alternatively 
to one of the criteria for the neutron skin proposed in 
Ref. [3l[ we consider another one which treats proton 
and neutron densities in a similar way. We extend the 
analysis of nuclear sizes presented in Ref. [22| by per- 
forming a more systematic study of a larger set of exotic 
nuclei and calculating also neutron skin thicknesses. The 
calculated charge rms radii are compared with the laser 
or muonic atoms spectroscopy measurements of isotope 
shifts performed on Sn [IE [II Sill, Ni El EH) and Kr 
[HI isotopes. The neutron skin thicknesses obtained in 
this paper are compared with the available experimental 
data extracted from methods mentioned above for even- 
even Sn isotopes with masses from 112 to 124. We also 
study whether the emergence of a skin is influenced by 



the nuclear shape, an issue that has not been sufficiently 
studied so far. The question of skin formation in nuclei 
having a non-spherical shape is discussed in detail on the 
example of Kr isotopes, assuming axial symmetry. 

The paper is organized in the following way. Section II 
contains the formalism of the deformed Skyrme HF+BCS 
method that provides the model density distributions, 
form factors, and nuclear radii. The numerical results 
and discussions are presented in Sec. III. Finally, we 
draw the main conclusions of this study in Sec. IV. 



II. DEFORMED SKYRME HF+BCS 
FORMALISM 

The results discussed in the next sections have been ob- 
tained from self-consistent deformed Hartree-Fock calcu- 
lations with density dependent Skyrme interactions (44| 
and pairing correlations. Pairing between like nucleons 
has been included by solving the BCS equations at each 
iteration either with a fixed pairing gap parameter (deter- 
mined from the odd-even experimental mass differences) 
or with a fixed pairing strength parameter. 

We consider in this paper the Skyrme force SLy4 [45| . 
We also show in some instances the results obtained from 
other parametrizations, namely Sk3 [4(| and SG2 (4?| be- 
cause they are among the most extensively used Skyrme 
forces and are considered as standard references. 

Assuming time reversal, the single-particle Hartree- 
Fock solutions for axially symmetric deformed nuclei are 
characterized by the eigenvalue fli of the third compo- 
nent of the total angular momentum on the symmetry 
axis and by the parity TTj. The state i can be written as 



$i (R,a,q) = XgM 



$+(r,z)e* A ^ X+ M 



+ <F-(r,z)e* A+ ^-M 



(1) 



where x 9i ((z), X±( a ) are isospin and spin functions, = 
fii ± 1/2 > 0. r, z, (p are the cylindrical coordinates of R. 

The wave functions "3^ are expanded into the 
eigenf unctions, 4> a , of an axially symmetric deformed 
harmonic-oscillator potential in cylindrical coordinates. 
We use 12 major shells in this expansion, 



(2) 



$i [R, cr, qj = x Qi (?) C ^ a °j ' 
with a — {n r , n z , A, £} and 

( - \ A eiA<P 

<j) a [R,aj = Tp nt Xr)^ nz (z)-j= Xs (v) , 
in terms of Hermite and Laguerre polynomials 
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<.(r) = 



(n r + A) 



/3 x ^^/2 e -V2 L A (??)) (5) 



with 



A = {mujz/h) 1 ' 2 , /3 X = (mwi/S) 1 / 2 , 

£ = , 7 7 = r 2 /3i. (6) 



A. Density distributions and root mean square 
radii 



The spin-independent proton and neutron densities are 
given by 



p(R) = p(r, z)=^2 2v iPi(r, z ) 



(7) 



in terms of the occupation probabilities v\ resulting from 
the BCS equations and the single-particle densities pi 



Pi (R) = Mr, z) = |*+(r, z)| 2 + \$i(r, z)\ 2 , (8) 



with 



£ ^,±1/2 ^ W (*) • (9) 



The multipole decomposition of the density can be writ- 
ten as [13, Eli 



p(r,z) = Ep A (i?)F A (cos0) 

A 

with multipole components A 



(10) 



px(R) = 



2A + 1 



Pa (cos t 



x p(R cos 9, Ram9)d(cos( 
and normalization given by 



(11) 



p(R)dR = X; 4tt / R dRpo(R) — X , (12) 



with X — Z, N for protons and neutrons, respectively. 

The mean square radii for protons and neutrons are 
defined as 



2 fR 2 Pp , n (R)dR 

^ ' p,n ^ 



/ p p , n (R)dR 



(13) 



and the rms radii for protons and neutrons are simply 
given by 



r =< r 2 > 1/2 

' P,n — ^ ' p,n ^ 



(14) 



The mean square radius of the charge distribution in 
a nucleus can be expressed as 



< r 2 > + < r c 2 h > p +(N/Z) < r 2 h > n 



+ r CM + r SO 



(15) 



where < r 2 > is the mean square radius of the point 
proton distribution in the nucleus (fl"5)) . < r 2 h > p and 
< r 2 h > n are the mean square charge radii of the charge 
distributions in a proton and a neutron, respectively. 
r^ M is a small correction due to the center of mass 
motion, which is evaluated assuming harmonic-oscillator 
wave functions. The last term rg G is a tiny spin-orbit 
contribution to the charge density. Correspondingly, we 
define the charge rms radius 



=< d > 1/2 



(16) 



B. Form factors and diffraction parameters 

Besides the mean square radii, additional characteris- 
tics of the density distributions can be deduced from the 
Fourier transforms of these densities. The form factors 
are defined as 



f p Aq) 



J p^ n (R)e^ R dR 
J p p , n (R)dR 



(17) 



In the Plane Wave Born Approximation (PWBA) the 
elastic electron scattering cross sections are related to the 
Fourier transform of the charge density 



F cb {q) = ^ I p ch {R)e^ A dR. 



(18) 



where q is the momentum transfer by the virtual photon 
in the scattering process. 

For each density multipole A, one defines a CA form 
factor 



A />00 

FCX (l) = Yj R^Rp x {R)j x {qR) . (19) 
In particular, 

F C0 (q) = YJ Q R 2 dRpo(R)jo(qR) (20) 



has the limit at q — > 

F C0 (q) - y/ R 2 dRpo(R) = l. (21) 



Elastic and inelastic electron scattering have been ex- 
tensively used to extract the various multipoles of the 
charge density, which show up in different transitions. In 
particular, in even-even deformed nuclei, F co (and corre- 
spondingly po) show up in the elastic cross section, while 
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F C2 (and hence p%) show up in the inelastic cross section 
for the transition + — > 2 + between the band-head and 
first excited rotational state [49l . I50I ] . 

In the next sections we will study the neutron skin 
thickness. We will use first the difference between the 
neutron and proton rms radii to characterize the different 
spatial extensions of neutron and proton densities. But 
as already noticed [201 ] . the rms radii (second moments 
of the densities) provide a very limited description of the 
nucleon density distributions. A more effective tool to 
analyze skins [2(| [Hj] is the Helm model [HI, [33| . This is 
a model that allows one to extract from the form factor in 
a simple way the two main characteristics of the density, 
a diffraction radius and a surface thickness. In this model 
one describes the density by convoluting a hard sphere 
(hs) density having diffraction radius Rd with a gaussian 
of variance cr, 



Taking out the factors which relate the rms radii 

to the radii of the equivalent uniform hard spheres, we 
define 

Rhs = < r 2 > 1 / 2 (29) 

and 

i?Hclm - ^573 i?™m = ^/i?2 +5(T 2 . (30) 

From these definitions we construct the following 
neutron-proton radius differences that will be used in the 
next sections 

AR d = R d (n) - R d (p) , (31) 



PHeim(r; Rd, cr) = phs(r; Rd) * Pc{r] cr) , (22) 



where 



and 



Phs(r,R d ) 



3X 
4ttR 3 „ 



e{R d -r), 



(23) 



p G (r;a) = (2^ 2 )- 3 / 2 e(-'- 2 / 2CT2 ). (24) 
The corresponding Helm form factor is 

Fndm{q) = F hs (q;R d )F G (q;a) 

= -^3i(qRd)e-° 2q2/2 . (25) 

qRd 

Now, the most prominent feature of the density distri- 
bution, namely its extension, can be related to the first 
zero in the form factor, this is the diffraction radius 



R d = 4.49341/gi 



(26) 



where q\ is the first zero of the form factor. 

The nuclear surface width a can be related to the 
height of the second maximum of the form factor located 

&f Qmax ■ 



■hi 



'In 



3jl(gmax-Rd) 
RdQ max 

F(q 

max 



(27) 



The variance cr is related to the surface thickness t (de- 
fined as the distance over which the density decreases 
from 90% to 10% of the central value) by t = 2.54 cr. 
Moreover, the surface thickness t is also related to the 
diffuseness a in the two-parameter Fermi distribution, by 
t = 4aln3 = 4.39 a. 

Taking into account that the second moment of a con- 
voluted distribution is given by the sum of the second mo- 
ments of the two single distributions, one gets the Helm 
rms radius 



R 



Hell! 



(i? 2 + 5o*) 



(28) 



AR hs = R hs (n 



- Rh S (p) 



< r 2 > x / 2 



(32) 



(33) 



Ai?Hclm — ^Holm(^) — ^Holm(p) • 



III. RESULTS AND DISCUSSION 



A. Root mean square radii and density 
distributions 



We start by showing our results for the rms radii of the 
charge distributions [Eq. [)16p]. We compare them to the 
available experimental information obtained from various 
methods including laser and muonic atoms spectroscopy 
[13, H [U HJ El |H, SI- We also compare our re- 
sults with different theoretical calculations. They include 
RMF calculations with NL3 parametrization and pairing 
correlations in BCS approach (RMF in Fig. [T|) [26j, 11011- 
relativistic calculations performed within HFB approach 
deduced under triaxial symmetry from D1S Gogny effec- 
tive interaction (HFB in Fig. [T]) , as well as calculations 
performed within a configuration mixing approach in the 
space spanned by the constrained HFB states. The lat- 
ter are done within the GCM under Gaussian overlap ap- 
proximation for the complete quadrupole collective space 
(GCM in Fig. CO) [13. 

Beginning with Sn isotopes for which more data and 
calculations are available, we show on the right panel 
of Fig. [T] our results for the squared charge radii differ- 
ences in Sn isotopes obtained from three different Skyrme 
forces, SLy4, SG2 and Sk3. We compare them to exper- 
iment, taking the radius of 120 Sn as the reference [4(j |. 
On the left panel we compare our SLy4 results for the 
charge radii with the other theoretical approaches men- 
tioned above. The general purpose of Fig. Q] is firstly 
to show that different Skyrme forces do not differ much 
in their predictions of charge rms radii and secondly, to 
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show that our results with SLy4 are comparable to other 
theoretical predictions including approaches that go be- 
yond the mean-field approximation, as well as relativis- 
tic approaches. Then, by comparing our results with 
experiment and with other theoretical results, we have 
evaluated the quality of our calculations. We conclude 
that our method reproduces the experimental data with 
a similar accuracy to other microscopic calculations that, 
as explained above, may be more sophisticated but may 
also be more time consuming. This agreement provides a 
good starting point to make predictions for other quan- 
tities such as neutron-proton radii differences, where the 
experimental information is scarce and it is not as accu- 
rate as in the case of charge radii. 

We complete this comparison of charge radii in Fig. [2] 
On the left panel we show our results for Ni isotopes and 
compare them with experiment (4ll . |42| and with results 
from RMF calculations [26] . On the right panel we show 
the same comparison for Kr isotopes. Data are taken 
from [43|| . In the Ni isotopes, we can see that the lower 
values of the rms radii occur around the double magic nu- 
cleus N = Z = 28, and around the semi-magic N = 50 in 
Kr isotopes. It is also worth mentioning that the bump 
shown around A — 76 in the RMF calculations of Kr iso- 
topes has its origin in the change of the ground-state nu- 
clear shape from oblate to prolate. In our case we obtain 
a smooth line because we only consider oblate shapes in 
this figure, as they correspond to the equilibrium shapes 
in most cases. 

Once we have confirmed that the agreement between 
our calculations with the experimental r c radii is satisfac- 
tory, we have guarantees that meaningful results will be 
obtained for the neutron and proton mean square radii 
(|14[) by using the same formalism with the same forces. 
Figure[3]contains our results with the SLy4 force for those 
radii in the three isotopic chains. They are compared 
with the predictions from RMF [26;]. We see that the 
tendency in the radii as a function of the mass number 
A is quite similar in both approaches, but in general the 
proton rms radii with Skyrme are systematically larger 
than the results from RMF. The situation is the oppo- 
site with respect to the neutron rms radii. At the same 
time the latter increase more slowly when calculated with 
SLy4. As a result we will get systematically differences 
between the neutron and proton rms radii, which are 
larger in the case of RMF as compared to the case of 
Skyrme forces. This is clearly seen in Fig.|4]where we plot 
the differences between the rms of neutrons and protons 
Ar np = r n — r p . On the left panel we show our results 
for Sn isotopes and compare them to RMF results and 
to experimental data taken from (p, p) scattering Q , 
antiprotonic atoms [icj |. giant dipole resonance method 
[IH , and spin dipole resonance method [HI, [l6| ■ As we 
can see in Fig. 2] the experimental data are located be- 
tween the predictions of both theoretical approaches and 
in general, there is agreement with experiment within the 
error bars. On the right panels we see the predictions for 
Ar np in the cases of Ni and Kr isotopes, where there are 



no data. 

The RMF results for the difference Ar np systematically 
overestimate the Skyrme HF results, as it can be seen 
from Fig.|U The reason for this is related to the difference 
in the nuclear symmetry energy and, consequently, to the 
different neutron equation of state (EOS) which has been 
extensively studied in recent years [5l|, [52, HH, HH . It was 
shown that there exists a linear correlation between the 
derivative of the neutron EOS (or the pressure of neu- 
tron matter) and the neutron skin thickness in heavy 
nuclei (defined as Arvp = r n — r p ) in both Skyrme HF 
[H H(§ and RMF [Hfl H3 models. We note that also 
a relation between Ar np and both volume and surface 
symmetry energy parameters was established recently by 
Danielewicz [58] and Steiner et al. [5jl] which provides 
a consistent description of nuclei with neutron excess. 
Typel and Brown j57[ demonstrated that the relativistic 
models produce larger neutron radii compared with the 
nonrelativistic ones, reflecting the fact that the satura- 
tion density of asymmetric matter is lower in the EOS 
when phenomcnological nucleon interaction in the RMF 
theory is used [60j. The results shown for neutron radii 
in Fig. 3 and correspondingly for neutron thicknesses in 
Fig. 4 support the above general conclusion. 

In the next figures we show the proton and neutron 
density distributions po(R) dTUJ) of some selected iso- 
topes in the three chains considered. We have chosen 
two extreme neutron-deficient and neutron-rich isotopes 
and one stable isotope between them. Figure [5] shows 
the neutron (solid) and proton (dashed) densities in the 
ioo,i20,i36g n isotopes. From left to right we see the evolu- 
tion of these densities as we increase the number of neu- 
trons. In the case of 100 Sn (iV=Z=50) we see that the 
two densities are practically the same except for Coulomb 
effects that make the protons to be more extended and, 
therefore, this has to be compensated with a small de- 
pression in the interior. The effect of adding more and 
more neutrons is to populate and extend the neutron den- 
sities. This makes also the proton distribution to follow 
the neutron one, increasing its spatial extension. The 
cost of this radius enlargement in the case of protons is a 
depression in the nuclear interior to preserve the normal- 
ization to the constant number of protons Z = 50. Then, 
it can be seen graphically the emergence of a region at the 
surface where the protons have practically disappeared 
while the neutrons still survive. We will quantify later 
this region in terms of the neutron skin thickness defini- 
tions. Figures [5] and [7] show the same information as in 
Fig. [5] but for 50,64js Ni and 70,84,98 Kr isotopes, respec- 
tively. The behavior of these densities corroborates the 
comments made on the case of Sn isotopes. 

As we mentioned in the last chapter, we will also char- 
acterize the skin thickness in terms of diffraction param- 
eters Rd and a deduced from the form factors. Figure [5] 
contains these form factors (fl7|) for protons and neutrons 
of the three Sn isotopes shown in Fig. [5] We can see how 
the diffraction zeroes at q\ (|26p and the location and 
magnitude of the second maximum, (7 max and F(q max ) 
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FIG. 1: Le ft p anel: Charge rms radii r c of tin isotopes. The SLy4 result is compared with the results from RMF calculations 
[H, HFB [|| and GCM [||. Experimental data are from [37|, H, IS S3 5 Ri S ht panel: Theoretical (with different Skyrme 
forces) and experimental isotope shifts 8{r%) of tin isotopes relative to 120 Sn. 
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(entering in Eq. (|27|) ) needed to extract Rd and the sur- 
face width a, change with the neutron number. Thus, 
we see that the q\ values diminish with increasing neu- 
tron number and therefore Rd increases accordingly for 
both protons and neutrons. The values of g max are also 
reduced when A increases but the values of the form fac- 
tor at these g max are rather similar. Consistently, the 
parameters a extracted from (|27j) are fairly similar. 



B. Neutron skin thickness 

The thickness of a neutron skin in nuclei may be de- 
fined in different ways. One of these possibilities is to 
define it as the difference between the root mean square 



radius of neutrons and that of protons, as we have plot- 
ted in Fig. 31 Similarly, it can be defined as the differ- 
ence between the neutron and proton radii of the equiv- 
alent uniform spheres [Eq. (|32pj. Alternatively, it can 
be defined as the difference between the neutron and 
proton diffraction radii (pTT| or Helm radii ([55)) . All of 
these quantities have already been discussed and used in 
the past as possible ways to quantify the skin thickness 
(see for example |20|), arriving to the conclusion that the 
radii difference defined in (|32p contains contribution from 
halo effects and the radii difference defined in (|33)) is a 
better measure of the skin. Nevertheless, qualitatively 
the difference between the two definitions becomes only 
apparent when dealing with very neutron- rich isotopes, 
which are presently beyond the experimentally observed 




isotopes and out of the scope of this paper. 

On the other hand, the skin thickness can be also de- 
fined in terms of some criteria that the neutron and pro- 
ton densities must fulfill. In Ref. [3l[ the neutron skin 
thickness is defined as the difference between two radii, 
Ri and R 2 . R\ is the radius at which the ratio of the 
neutron density to the proton density is equal to some 



given value (4 in [3l[). R? is the radius at which the 
neutron density becomes smaller than some perce ntag e 
of the density at the center of the nucleus (1 % in [3l|). 
When this difference, AR = R2 — Ri, is larger than some 
established value (in [3l| this value is 1 fm, which is com- 
parable to the range of the nuclear force) , a neutron skin 
with skin thickness AR is said to occur. 
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FIG. 4: Difference between neutron and proton rms radii Ar np of Sn, Ni, and Kr isotopes calculated with SLy4 force. The 
RMF calculation results are from Ref. [2q|. The experimental data for Sn isotopes measured in (p,p) reaction (open stars) 
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FIG. 5: HF+BCS proton and neutron densities p (R) of 100 Sn, 120 Sn, and 136 Sn calculated with SLy4 force 
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FIG. 8: Proton (left panel) and neutron (right panel) form factors for the 100 Sn, 120 Sn, and 136 Sn isotopes are calculated in 
the PWBA. 



The factors used to define the skin thickness in the 
above criteria could have been differently chosen in rather 
arbitrary manners. Therefore, the absolute sizes of the 
skin thickness do not have a very precise meaning. Nev- 
ertheless, these values are useful to judge how the nucleon 
skins develop as the number of nucleons change. Indeed, 
we have also considered the case where the first criterion 
for the inner radius R\ of the neutron skin is changed. 
We use instead of the above criterion for the radius 
at which the proton density becomes smaller than 1% of 
the latter at the center, which is similar to the criterion 
used to define the outer radius i?2, but in this case for 
proton density instead of the neutron density. When we 
use the conditions in Ref. [3l| . we call it criterion (a). 
When we use the alternative condition for R\ , we call it 
criterion (b). 

We show in Fig. [5] the results obtained for the neu- 
tron skin thickness in Sn isotopes according to the differ- 
ent definitions discussed above. The left panel contains 
the results for definitions involving directly the difference 
between neutron and proton radii, either the equivalent 



hard spheres radii Ai?h s [Eq. (|32)) ] corresponding to the 
rms radii, the diffraction radii Ai?d [Eq. J5T|) ] . and the 
Helm radii Ai?Hcim [Eq. (JS31)] • The skin thickness pre- 
dicted by the difference of the very simple diffraction radii 
is in general smaller than the thickness predicted by the 
other two more involved options that are very similar in 
this range of masses. The right panel contains the neu- 
tron skin thickness defined according to the criteria on 
the density distributions (a) (solid line) and (b) (dashed 
line). They only differ in the way in which the starting 
radius of the skin R\ is chosen. One can see that we ob- 
tain larger neutron skin thicknesses when using criterion 
(b) in the lighter isotopes, but this is reversed for heavier 
isotopes and we get larger thickness when using criterion 
(a). This fact is confirmed also by the values of the radii 
Ri and R2 and their differences AR listed in Table [T] for 
the heaviest three isotopes in each chain considered. In 
general, the formation of a skin when using (a) starts at 
distances smaller than those in case (b) or comparable 
with them, which leads to larger absolute size of the neu- 
tron skin produced by criterion (a). It is in this region of 
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heavier isotopes where we can properly talk about a neu- 
tron skin formation. In this region, criterion (b) somehow 
establishes a lower limit for the skin thickness. The latter 
can be arbitrarily enlarged by relaxing the p n /p p condi- 
tion to values lower than 4. Similar comments apply also 
to the next figures, Fig. [10] for Ni isotopes and Fig. [ITJ 
for Kr isotopes. 

TABLE I: Radii R[ a) , R{ b) , R 2 and their differences (skin 
thicknesses) Ai? (a) = R 2 - R[ a) and AR W = R 2 - R[ b) (in 
fm) according to criteria (a) and (b). 



Nuclei 


R[ a) 




Ra 


AR (a) 


ARM 


V4 Ni 


6.49 


6.48 


7.20 


0.71 


0.72 


76 Ni 


6.38 


6.49 


7.24 


0.86 


0.75 


78 Ni 


6.32 


6.52 


7.29 


0.97 


0.77 


96 Kr 


7.32 


7.31 


7.83 


0.51 


0.52 


98 Kr 


7.20 


7.33 


7.90 


0.70 


0.57 


ioo Kr 


7.06 


7.34 


7.96 


0.90 


0.62 


132 Sn 


7.64 


7.76 


8.08 


0.44 


0.32 


134 Sn 


7.50 


7.78 


8.19 


0.69 


0.41 


136 Sn 


7.40 


7.82 


8.32 


0.92 


0.50 



We would like to emphasize that although different def- 
initions of the neutron skin thickness produce different 
absolute values for it, the relative skin thicknesses cor- 
responding to the evolution as the number of neutrons 
increase indicates the formation of such a skin that can 
be expected to start at A > 132 in Sn, A > 74 in Ni, and 
A > 96 in Kr isotopes, as it is observed in Figs. [9lfTTl 

Finally, we also consider the most neutron-deficient re- 
gion of Ni isotopes in a search for the formation of a 
proton skin. We have already seen in the left panel in 
Fig. [IHl that the neutron skin thickness defined in terms 
of differences between neutron and proton radii become 
negative at some point, indicating that the proton dis- 
tribution extends beyond the neutron one. This can be 
further explored by reversing the definitions of R\ and 
i?2 interchanging the role of protons and neutrons. We 
show the results in the inset of the right panel in Fig. [T0l 
where we have applied the criterion (b) with protons and 
neutrons interchanged. We find no proton skin when ap- 
plying criterion (a). One can see that a small skin starts 
developing in these isotopes but we cannot push it fur- 
ther because 48 Ni is already at the proton drip line. The 
results are then not conclusive enough to assess the exis- 
tence of a proton skin in these isotopes. This possibility 
could be explored in the future in the most proton-rich 
nuclei approaching the proton drip lines of lighter nuclei 
with Z > N. 



C. Neutron skin and deformation 

When the nucleus is deformed, the thickness of the 
neutron skin might depend on the direction. It is an 
interesting and natural question to ask whether the de- 
formed densities give rise to a different skin size in the 



different directions. It is also interesting to know whether 
the emergence of the skin may be influenced by the nu- 
clear shape. We study in this work such a dependence 
on the example of Kr isotopes, which are examples of 
well deformed nuclei characterized by a large variety of 
competin g nu clear shapes [6l|. Constraint HF+BCS cal- 
culations [6ll. l62l | show also the possibility of shape co- 
existence in these nuclei. The results which we obtain 
for the binding energy of the three previously selected 
Kr isotopes as a function of the quadrupole parameter 
f3 = n/bQp/ (Zrp) (Q p being the proton quadrupole 
moment) are presented in Fig. 1121 In this figure the dis- 
tance between two ticks in the vertical axis is always 1 
MeV but the origin is different for each curve. As we can 
see, both prolate and oblate shapes produce minima very 
close in energy. Then, we have chosen the neutron rich 
isotope 98 Kr to study the sensitivity of the neutron skin 
thickness to the various directions in the two shapes. 

GCM calculations built on the constrained HF+BCS 
states may be carried out in order to describe more prop- 
erly some ground-state properties in deformed nuclei. In 
the case of 98 Kr the potential energy curve (Fig. [T2"|) 
shows pronounced minima at oblate and prolate shapes, 
which are separated by an energy barrier of about 6 MeV. 
Thus, one expects the ground state of 98 Kr to be basically 
described by a linear combination of these two configu- 
rations. 

We first study the intrinsic density distributions p{R) 
in various selected directions. For that purpose we show 
in Figs. [13] and [14] the densities of 98 Kr for oblate and 
prolate shapes, respectively. We can see the spatial dis- 
tributions for neutrons (solid) and protons (dotted) in 
three different directions: z-direction (r = 0), r-direction 
(z = 0), and r = z direction. We can observe that the 
profiles of the densities as well as the spatial extensions 
change with the direction. Clearly, the densities are more 
extended in the z-direction in the case of prolate shapes. 
The opposite is true in the case of oblate shapes. The 
case r = z gives always intermediate densities. We have 
added in the three directions a couple of full dots, indi- 
cating the radii R\ and i?2 that defines the skin thickness 
according to the above mentioned criterion (a). 

The dependence of the intrinsic density on the different 
directions can be also seen in Fig. 1151 where we plot as an 
example the proton densities in the three directions men- 
tioned above for oblate (left) and prolate (right) shapes 
in the same plane. We see more clearly how the exten- 
sion of the density in the ^-direction (labeled r = 0) is 
the largest for the prolate shape and the shortest for the 
oblate shape. We also plot for comparison the monopole 
component po (R) (fTU)) that lies between the two extreme 
cases and it is close to the density in the r — z direction. 

It is also worth looking at the points in the (r, z) plane 
that define the ellipses where the criteria for R\ and Ri 
are met. Figure [TBI shows these points for protons (thin 
lines) and neutrons (thick lines) and for the two shapes, 
prolate (solid) and oblate (dashed). We can see that the 
size of the skin changes little with the directions perpen- 
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FIG. 11: Same as in Fig. [9] but for Kr isotopes. 
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FIG. 12: Binding energies E calculated with the SLy4 force as a function of the quadrupole parameter j3 for the even-even 
70 Kr, 84 Kr, and 98 Kr isotopes. 




FIG. 14: Same as in Fig. [13] but for prolate shape of 98 Kr. 
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dicular to the surface, but shows a tendency to increase 
on the shorter axis. It is interesting to note that the 
skin size of the spherical component po(R) is an inter- 
mediate value. The overall skin thickness is also simi- 
lar in the oblate and prolate equilibrium shapes. From 
this example we may conclude that the skin thickness 
does not depend much on the oblate or prolate character 
of the deformation. This is in line with the conclusions 
reached in Ref. [35| on the example of Dy isotopes, where 
it was shown that the neutron skin is nearly indepen- 
dent of the size of deformation (spherical, deformed or 
superdeformed) . 

Figure [T71 shows the monopole, po(R), and quadrupole, 

P2{R), components of the intrinsic density p(R) (JTOj) for 
protons (dashed lines) and neutrons (solid lines) and for 
the oblate and prolate shapes in 98 Kr. We can see that 
P2(R) is peaked at the surface positively in the case of 
the prolate deformation and negatively in the case of 
the oblate one. This makes the total density in the in- 
direction to be incremented with respect to the po den- 
sity in the prolate case and to be decreased in the oblate 
one. The opposite is true with respect to the direc- 
tion perpendicular to the symmetry axis z. We can also 
see that the skin thickness derived from the po compo- 
nents is quite similar to the thickness derived from the 
quadrupole components p2- This explains the approxi- 
mately constant skin thickness observed in the different 
directions in Fig. 1161 



IV. CONCLUSIONS 

In this work we perform a theoretical analysis of nu- 
clear skins, exploring various definitions. For this pur- 
pose we examine three chains of Ni, Kr, and Sn isotopes 
which might be of particular interest in the future ex- 
periments in GSI and RIKEN. The densities of these 



nuclei are calculated within a deformed HF+BCS ap- 
proach with Skyrme-type density-dependent effective in- 
teractions [6l|. We have shown that this model gives 
a very reasonable description of the charge rms radii of 
the Sn, Ni, and Kr isotopes and of the differences be- 
tween neutron and proton rms radii for several Sn iso- 
topes. This is confirmed by the good agreement with the 
available experimental data, as well as with other theo- 
retical predictions. Three Skyrme parametrizations have 
been involved in the calculations: SG2, Sk3 and SLy4. 
Most of the results shown in the paper are obtained with 
SLy4 force, but the other Skyrme interactions produce 
similar results. 

For the first time the various definitions which have 
been previously proposed to determine the neutron skin 
thickness, involving both matter radii and tails of nu- 
clear densities, have been compared within a deformed 
Skyrme HF+BCS model. We find that all definitions of 
the neutron skin predict to a different extent the exis- 
tence of a skin in nuclei far from the stability line. Par- 
ticularly, a pronounced neutron skin can be attributed to 
heavier isotopes of the three chains considered, namely 
with A > 132 for Sn, A > 74 for Ni, and A > 96 
for Kr isotopes. We also find that for a given isotopic 
chain the increase of the skin with the neutron number in 
the neutron-rich nuclei exhibits a rather constant slope, 
which is different depending on the definition of nuclear 
skin. More significant neutron skin is obtained when ana- 
lyzing its formation by means of definition from Ref. (3l| 
(called criterion (a)) or using an alternative one (called 
criterion (b)). In this case we get an absolute size of 
the skin larger than 0.4 fm and almost reaching 1 fm for 
the heaviest isotopes (in the case of criterion (a)). At 
the same time, the neutron skin determined by the dif- 
ference between neutron and proton radii using diffrac- 
tion parameters defined in the Helm model shows a more 
smooth gradual increase with the neutron excess and it is 
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FIG. 16: Radii _Ri and R2 according to criterion (a) for neutrons (thick lines) and protons (thin lines) in 98 Kr nucleus (shown 
in rz plane) corresponding to its oblate (solid lines) and prolate (dashed lines) shape. 




R [fm] R [fm] 

FIG. 17: Monopole po(R) (thin lines) and quadrupole p2{R) (thick lines) neutron (solid lines) and proton (dashed lines) density 
distributions of 98 Kr with oblate and prolate shape. 



in size of around 0.3-0.4 fm. We would like to note that 
our results for Sn isotopes are consistent with the results 
of calculations from Ref. [20( with SLy4 parametrization. 
In both calculations the analysis of neutron skin forma- 
tion is based on the nuclear form factors, which are well 
suited for such study since the diffraction parameters are 
mainly sensitive to the nuclear densities in the surface 
region. 

We also show on the example of the neutron-deficient 
Ni isotopes the possibility to find a proton skin in a simi- 
lar way to the neutron skin. Although the analysis, which 
was performed in our paper for this case, uses an alter- 
native criterion to that applied in [31( , it indicates a sit- 
uation close to proton skin formation in Ni isotopes very 



close to the proton drip line. However, the search for the 
existence of proton skin could be explored in the most 
proton-rich nuclei approaching the proton drip lines of 
lighter nuclei, where Z > N. 

In the present work the effects of deformation on the 
skin formation are studied in Kr isotopes which are well 
deformed nuclei. Taking as an example 98 Kr isotope, we 
find that the profiles of the proton and neutron densi- 
ties, as well as the spatial extensions change with the 
direction in both oblate and prolate shapes. At the same 
time, the neutron skin thickness remains almost equal 
along the different directions perpendicular to the sur- 
face. Same type of calculations have been also performed 
on the example of 100 Kr, exhibiting a similar potential en- 
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ergy curve. In this case, the conclusion concerning neu- 
tron skin thickness on the different directions remains 
unchanged. We find a very weak dependence of the neu- 
tron skin formation on the character of deformation. This 
is a useful information, worth to be known before com- 
plete GCM calculations are performed because it indi- 
cates that no drastic changes in the neutron skin thick- 
ness are expected when such more sophisticated calcula- 
tions are performed. 

The results obtained in the present paper demonstrate 
the ability of our microscopic theoretical method to pre- 
dict the nuclear skin in exotic nuclei. They also illustrate 
the range of the skin sizes to be expected depending on 
the adopted skin definition. More definite conclusions on 
the emergence of nuclear skin will be drawn when direct 
measurements of proton and neutron form factors, and 



thus the corresponding proton and neutron densities, for 
these nuclei will become available at the upcoming ex- 
perimental facilities. 
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